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A  study  was  performed  to  elucidate  the  mechanisms  of  extinction  and  autoignition  of  n -heptane  in 
strained  laminar  flows  under  non-premixed  conditions.  A  previously  developed  detailed  mechanism  made 
up  of  2540  reversible  elementary  reactions  among  556  species  was  the  starting  point  for  the  study.  The 
detailed  mechanism  was  previously  used  to  calculate  ignition  delay  times  in  homogeneous  reactors,  and 
concentration  histories  of  a  number  of  species  in  plug-flow  and  jet-stirred  reactors.  An  intermediate  mech¬ 
anism  made  up  of  1282  reversible  elementary  reactions  among  282  species  and  a  short  mechanism  made 
up  of  770  reversible  elementary  reactions  among  159  species  were  assembled  from  this  detailed  mecha¬ 
nism.  Ignition  delay  times  in  an  isochoric  homogeneous  reactor  calculated  using  the  intermediate  and  the 
short  mechanism  were  found  to  agree  well  with  those  calculated  using  the  detailed  mechanism.  The 
intermediate  and  the  short  mechanism  were  used  to  calculate  extinction  and  autoignition  of  n-heptane  in 
strained  laminar  flows.  Steady  laminar  flow  of  two  counterflowing  streams  toward  a  stagnation  plane  was 
considered.  One  stream,  made  up  of  prevaporized  n-heptane  and  nitrogen,  was  injected  from  the  fuel 
boundary,  and  the  other  stream,  made  up  of  air  and  nitrogen,  was  injected  from  the  oxidizer  boundary. 
Critical  conditions  of  extinction  and  autoignition  given  by  the  strain  rate,  temperature,  and  concentrations 
of  the  reactants  at  the  boundaries  were  calculated.  The  results  were  found  to  agree  well  with  experiments. 
Sensitivity  analysis  was  carried  out  to  evaluate  the  influence  of  various  elementary  reactions  on  autoignition. 
At  all  values  of  the  strain  rate  investigated  here,  high-temperature  chemical  processes  were  found  to  control 
autoignition.  In  general,  the  influence  of  low-temperature  chemistry  was  found  to  increase  with  decreasing 
strain.  A  key  finding  of  the  present  study  is  that  strain  has  more  influence  on  low-temperature  chemistry 
than  the  temperature  of  the  reactants. 


Introduction 

Chemical-kinetic  mechanisms  that  describe  oxi¬ 
dation  of  n-heptane  (n-C7H16)  have  been  developed 
by  numerous  investigators  [1-13].  n -Heptane  is  a 
primary  reference  fuel  for  octane  rating  in  internal 
combustion  engines.  It  has  a  cetane  number  of  ap¬ 
proximately  56,  which  is  similar  to  the  cetane  num¬ 
ber  of  conventional  diesel.  Recently,  a  comprehen¬ 
sive  chemical-kinetic  mechanism  was  developed  to 
describe  combustion  of  n-heptane  [10],  Numerical 
calculations  were  performed  using  this  mechanism, 
and  the  results  were  compared  with  experimental 
data  obtained  in  flow  reactors,  shock  tubes,  and 
rapid  compression  machines.  Here,  chemical-kinetic 
mechanisms  are  assembled  from  this  comprehensive 
mechanism  and  are  used  to  describe  extinction  and 
autoignition  of  n-heptane  in  non-premixed  systems. 
The  results  are  compared  with  experiments  con¬ 
ducted  in  the  counterflow  configuration  [12-14], 

The  present  study  was  motivated  by  a  need  to  de¬ 
velop  chemical-kinetic  mechanisms  that  can  be  used 


in  computational  fluid  dynamic  codes  to  model  com¬ 
bustion  and  transport  processes  taking  place  in  die¬ 
sel  engines.  Combustion  processes  in  diesel  engines 
closely  resemble  non-premixed  systems.  Previous 
experimental  studies  on  autoignition  were  focused 
on  measuring  ignition  delay  times  in  shock  tubes 
[15]  and  rapid  compression  machines  [16-18]. 
These  studies  were  conducted  on  premixed  systems 
and  revealed  various  aspects  of  n-heptane  autoigni¬ 
tion  and  combustion.  The  present  study  addresses 
extinction  and  autoignition  of  n-heptane  in  counter¬ 
flow  non-premixed  systems  where  transport  pro¬ 
cesses  play  an  important  role.  Computational  re¬ 
sources  required  for  calculating  aspects  of 
autoignition  in  a  one-dimensional  counterflow  con¬ 
figuration  are  significantly  larger  than  those  required 
in  a  homogeneous  system.  Large  reaction  mecha¬ 
nisms,  such  as  those  for  n-heptane,  need  to  be  re¬ 
duced  to  make  the  numerical  problem  tractable  in 
one-dimensional  configurations. 

Lindstedt  and  Maurice  [5]  have  developed  a  de¬ 
tailed  chemical-kinetic  mechanism  for  describing 
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LAMINAR  DIFFUSION  FLAMES— Counterflow  Flames 
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Fig.  1.  Schematic  illustration  of  the  counterflow 
configuration. 


the  oxidation  of  n-licptane.  This  mechanism  has 
been  tested  by  comparing  numerical  results  with  ex¬ 
perimental  data  obtained  in  stirred  reactors,  coun¬ 
terflow  non-premixed  flames,  and  premixed  flames. 
This  was  one  of  the  first  attempts  to  develop  a  chem¬ 
ical-kinetic  mechanism  for  n-heptane  which  can  be 
applied  to  flames  where  transport  processes  are  im¬ 
portant,  and  to  stirred  reactors.  This  mechanism, 
however,  did  not  address  ignition  and  extinction  of 
non-premixed  flames.  Seiser  et  al.  [12]  calculated 
the  structure  and  critical  conditions  of  extinction  of 
non-premixed  n-heptane  flames  and  compared  the 
results  with  experimental  measurements.  This  study 
employed  a  chemical-kinetic  mechanism  which  had 
been  previously  simplified  by  introducing  ad  hoc  ap¬ 
proximations  [8].  At  given  concentrations  of  the  re¬ 
actants,  the  calculated  strain  rate  at  extinction  was 
found  to  be  higher  than  that  measured  [12],  Here, 
critical  conditions  of  extinction  are  calculated  using 
mechanisms  that  were  deduced  from  the  compre¬ 
hensive  mechanism  of  Curran  et  al.  [10],  and  the 
results  are  compared  with  experimental  data  shown 
in  Ref.  [12],  Experimental  data  for  the  critical  con¬ 
ditions  of  autoignition  of  n-heptane  in  the  counter¬ 
flow  configuration  are  given  in  Refs.  [13,14],  Using 
these  experimental  data,  overall  chemical-kinetic 
rate  parameters  that  characterize  the  rate  of  one- 
step  overall  reaction  between  fuel  and  oxygen  were 
obtained  [14],  Here,  critical  conditions  of  autoigni¬ 
tion  are  calculated  using  mechanisms  that  were  de¬ 
duced  from  the  comprehensive  mechanism  [  10] ,  and 
the  results  are  compared  with  experimental  data 
shown  in  Refs.  [13,14], 


the  counterflow  configuration  employed  in  previous 
experimental  studies  [12-14]  and  in  the  present  nu¬ 
merical  study.  In  this  configuration,  a  fuel  stream 
made  up  of  prevaporized  n-heptane  and  nitrogen  is 
injected  from  one  duct,  called  the  fuel  duct,  and  an 
oxidizer  stream  made  up  of  air  and  nitrogen  is  in¬ 
jected  from  the  other  duct,  called  the  oxidizer  duct. 
These  jets  flow  into  the  mixing  layer  between  the 
two  ducts.  The  exit  of  the  fuel  duct  is  the  fuel  bound¬ 
ary,  and  the  exit  of  the  oxidizer  duct  is  the  oxidizer 
boundary.  The  mass  fraction  of  fuel,  the  tempera¬ 
ture,  and  the  component  of  the  flow  velocity  normal 
to  the  stagnation  plane  at  the  fuel  boundary  are  rep¬ 
resented  by  Yf1,  T1 ,  and  V1;  respectively.  The  mass 
fraction  of  oxygen,  the  temperature,  and  the  com¬ 
ponent  of  the  flow  velocity  normal  to  the  stagnation 
plane  at  the  oxidizer  boundary  are  represented  by 
Yo2,2-  t2,  and  V2,  respectively.  The  tangential  com¬ 
ponents  of  the  flow  velocities  at  the  boundaries  are 
presumed  to  be  equal  to  zero.  The  distance  between 
the  fuel  boundary  and  the  oxidizer  boundary  is  rep¬ 
resented  by  L. 

In  the  experiments  and  numerical  calculations,  the 
momenta  of  the  counterflowing  reactant  streams 
PiVf,  i  =  1,  2  at  the  boundaries  are  kept  equal  to 
each  other.  Here,  /q  and  p2  represent  the  density  of 
the  mixture  at  the  fuel  boundary  and  the  oxidizer 
boundary,  respectively.  This  condition  ensures  that 
the  stagnation  plane  formed  by  the  two  streams  is 
approximately  in  the  middle  of  the  region  between 
the  two  boundaries.  The  value  of  the  strain  rate,  de¬ 
fined  as  the  normal  gradient  of  the  normal  compo¬ 
nent  of  the  flow  velocity,  changes  from  the  fuel 
boundary  to  the  oxidizer  boundary  [19].  The  char¬ 
acteristic  strain  rate  on  the  oxidizer  side  of  the  stag¬ 
nation  plane,  a2,  is  presumed  to  be  given  by  [19] 


(1) 


Equation  1  is  obtained  from  an  asymptotic  theory  in 
which  the  Reynolds  numbers  of  the  laminar  flow  at 
the  boundaries  are  presumed  to  be  large  [19],  Criti¬ 
cal  conditions  of  extinction  are  presumed  to  be  given 
by  the  strain  rate,  a2  e,  the  temperatures,  and  the 
mass  fraction  of  fuel  and  oxygen  at  the  boundaries. 
Critical  conditions  of  autoignition  are  presumed  to 
be  given  by  the  strain  rate,  a2>1,  the  temperature  of 
the  oxidizer  stream,  T2I,  the  temperature  of  the  fuel 
stream,  and  the  mass  fraction  of  fuel  and  oxygen  at 
the  boundaries.  The  experiments  were  conducted  at 
a  pressure  of  1.013  bar. 


Description  of  Experimental  and  Numerical 
Studies  in  the  Counterflow  Configuration 

Steady,  axisymmetric,  laminar  flow  of  two  coun- 
terflowing  streams  toward  a  stagnation  plane  was 
considered.  Fig.  1  shows  a  schematic  illustration  of 


Experimental  Procedure 

A  detailed  description  of  the  burners  has  been 
given  elsewhere  [12-14],  The  flow  rates  of  gases 
were  measured  by  computer-regulated  mass  flow 
controllers.  The  velocities  of  the  reactants  at  the 
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boundaries  were  presumed  to  be  equal  to  the  ratio 
of  their  volumetric  flowrates  to  the  cross-sectional 
area  of  the  ducts.  The  temperature  of  the  fuel  stream 
and  the  temperature  of  the  oxidizer  stream  at  the 
boundaries  were  measured  using  thermocouples.  A 
brief  description  of  the  experimental  procedure  is 
given  here. 

Critical  Conditions  of  Extinction 

Extinction  experiments  were  previously  carried 
out  allowing  only  small  changes  in  the  flame  position 
in  the  reactive  flow  field  [12],  It  is  convenient  to 
express  the  flame  position  in  terms  of  a  conserved 
scalar  quantity  £,  called  the  mixture  fraction.  The 
mixture  fraction  is  so  defined  that  =  1.0  in  the 
fuel  stream  and  Q  =  0  in  the  oxidizer  stream  [20], 
The  location  of  the  flamesheet,  <Jst,  where  the  flux  of 
the  fuel  and  the  flux  of  oxygen  are  in  stoichiometric 
proportion,  is  given  by  <Jst  =  [1  +  11Yf1Wq2/ 
(Y02>2Wf)]  ~~  \  where  WF  and  WQo  represent  the  mo¬ 
lecular  weights  of  fuel  and  oxygen,  respectively.  Ex¬ 
tinction  experiments  were  performed  keeping  fst 
constant  at  0.1.  The  temperature  of  the  fuel  stream 
was  Tf  =  345  K,  and  the  temperature  of  the  oxidizer 
stream  was  T2  =  298  K.  The  distance  between  the 
fuel  boundary  and  the  oxidizer  boundary  was  L  = 
10  mm.  At  some  selected  value  of  Yq22>  the  flame 
was  stabilized  at  «2  <  aZe.  The  strain  rate  was  in¬ 
creased  by  increasing  V1  and  V2  until  extinction  was 
observed.  Experimental  results  are  shown  later. 

Critical  Conditions  of  Autoignition 

Previous  autoignition  experiments  were  con¬ 
ducted  with  the  mass  fraction  of  prevaporized  fuel, 
Yf1,  maintained  at  0.387  [13,14],  The  temperature 
at  the  fuel  boundary,  Th  was  378  K.  The  oxidizer 
stream  was  air  with  a  mass  fraction  of  oxygen 
Yo2,2  =  0.233.  The  distance  between  the  fuel 
boundary  and  the  oxidizer  boundary  was  L  =  12 
mm.  At  a  given  strain  rate  and  oxidizer  temperature, 
T*2  <  T2  I,  the  flow  field  was  established.  The  tem¬ 
perature  at  the  oxidizer  boundary  was  gradually  in¬ 
creased  until  autoignition  took  place.  Experimental 
results  are  shown  later. 

Formulation  of  the  Numerical  Problem 

The  conservation  equations  of  mass,  momentum, 
and  energy  and  the  species  balance  equations  used 
in  the  formulation  of  the  numerical  problem  have 
been  summarized  elsewhere  [8,21,22],  The  species 
balance  equations  include  thermal  diffusion,  and  the 
energy  conservation  equation  includes  radiative  heat 
loss  from  carbon  dioxide  and  water  vapor  [22],  Buoy¬ 
ancy  is  neglected.  Calculations  are  performed  over 
a  computational  domain  of  10  mm  to  obtain  the  criti¬ 
cal  conditions  of  extinction  and  over  a  computational 


domain  of  12  mm  to  obtain  the  critical  conditions  of 
autoignition.  At  both  ends  of  the  computational  do¬ 
main,  the  mass  fractions  of  the  reactants  and  the 
normal  components  of  the  flow  velocity  are  speci¬ 
fied.  The  values  of  the  tangential  component  of  the 
flow  velocity  at  both  ends  are  set  equal  to  zero  (the 
so-called  plug-flow  boundary  conditions).  Critical 
conditions  of  extinction  are  calculated  with  Tl  =  345 
K  and  T2  =  298  K.  Critical  conditions  of  autoigni¬ 
tion  are  calculated  with  T1  =  378  K.  The  character¬ 
istic  strain  rate  at  the  stagnation  plane  is  calculated 
using  equation  1. 

The  chemical-kinetic  mechanism  used  to  calculate 
the  critical  conditions  of  extinction  and  autoignition 
was  deduced  from  the  comprehensive  mechanism  of 
Curran  et  al.  [10],  This  mechanism,  called  the  de¬ 
tailed  mechanism,  comprises  2540  reversible  ele¬ 
mentary  reactions  among  555  species.  The  rate  con¬ 
stant,  k,  of  any  elementary  reaction  in  this 
mechanism  is  given  by  k  =  BV 1  exp[  —  E/CRT)], 
where  T  is  the  temperature  and  R  is  the  gas  con¬ 
stant.  The  frequency  factor,  B,  the  temperature  ex¬ 
ponent,  n,  and  the  activation  energy,  E ,  are  the  rate 
parameters.  This  comprehensive  mechanism  was 
previously  used  to  calculate  ignition  delay  times  in 
homogeneous  reactors,  and  concentration  histories 
of  a  number  of  species  in  plug-flow  and  jet-stirred 
reactors.  The  results  were  found  to  agree  with  ex¬ 
periments  [10],  Improvements  have  been  made  to 
this  mechanism  [23],  Rate  parameters  for  the  for¬ 
ward  rate  constant  of  alkyl  radical  R  addition  to  ox¬ 
ygen  represented  by  the  elementary  reaction  R  + 
02  =  R02  is  now  given  by  B  =  4.52  X  1012  cm3/ 
(mol  s),  n  =  0,  andE  =  0  for  the  primary,  B  —  7.54 
X  1012  cm3/(mol  s),  n  =  0,  and  E  =  0  for  the 
secondary,  and  B  =  1.41  X  1013  cm3/(mol  s),  n  — 
0,  and  E  —  0  for  the  tertiary  alkyl  radical.  With  the 
alkyl  denoted  by  R',  rate  parameters  for  the  forward 
rate  constant  of  the  reaction  R'Oa  +  R  =  R'O  + 
RO  are  B  =  7.00  X  1012  cm3/(mol  s),  n  =  0,  and 
E  =  —4184  J/mol;  for  the  reaction  R02  +  H02  = 
R02H  +  02,  B  =  1.75  X  1010  cm3/(mol  s),  n  = 
0,  and  E  —  —  13702.58  J/mol;  and  for  the  reaction 
ROo  +  R'02  —  02  +  RO  +  R'O,  B  —  1.40  X 
1016  cm3/ (mol  s),  n  =  -  1.61,  and  E  =  7782.23  J/ 
mol.  Modifications  in  reaction  rate  rules  were  doc¬ 
umented  in  Ref.  [23],  Recent  changes  in  the  ther¬ 
modynamic  properties  of  peroxy  radicals  have  been 
adopted  [24,25], 

The  number  of  reactions  and  species  in  the  de¬ 
tailed  mechanism  was  reduced  before  it  was  used  to 
calculate  the  critical  conditions  of  extinction  and  au¬ 
toignition  in  non-premixed  systems.  To  reduce  the 
mechanism,  ignition  delay  times,  Tjg,  in  an  isochoric 
homogeneous  reactor  configuration  were  calculated 
using  the  detailed  mechanism.  Ten  conditions  were 
selected,  given  by  values  of  initial  pressure  equal  to 
1  bar  and  13.5  bar  and  values  of  initial  temperatures 
equal  to  625  K,  740  K,  909  K,  1176  K,  and  1667  K. 
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1000  /r  [1/K] 

Fig.  2.  Ignition  delay  time,  rig,  in  stoichiometric 
mixtures  of  n-heptane  vapor  and  air  as  a  function  of  the 
reciprocal  of  initial  temperature,  Tu,  for  various  pressures. 
The  short-dashed  line  represents  rig  calculated  using  the 
detailed  mechanism.  The  long-dashed  line  represents  rig 
calculated  using  the  intermediate  mechanism.  The  solid 
line  represents  rig  calculated  using  the  short  mechanism. 
The  symbols  represent  experimental  data  obtained  in  shock 
tubes  [15]  and  rapid  compression  machines  [16], 


Fig.  3.  The  maximum  temperature,  Tmax,  as  a  function 
of  the  strain  rate.  Calculations  are  performed  using  the 
short  mechanism  for  73  =  1155  K  and  =  378  K.  The 
mass  fraction  of  fuel  in  the  fuel  stream  is  YF1  =  0.387. 

To  test  the  influence  of  any  given  species  on  the 
ignition  delay  time,  all  elementary  reactions  in  which 
this  species  appears  were  removed  from  the  detailed 
mechanism,  and  the  values  of  rig  were  computed  at 
the  selected  10  conditions  using  the  truncated  de¬ 
tailed  mechanism  now  made  up  of  555  species.  The 
deviation  of  the  ignition  delay  time  calculated  using 
this  truncated  mechanism  from  the  ignition  delay 
time  calculated  using  the  detailed  mechanism  was 


obtained.  If  the  deviation  at  all  conditions  was  less 
than  0.84%,  this  species  was  removed  from  the 
mechanism.  Some  species  were  retained  in  the 
mechanism  even  if  the  deviation  was  less  than 
0.84%,  because  their  influence  on  critical  conditions 
of  extinction  and  autoignition  was  suspected  to  be 
greater  than  their  influence  on  rig  in  homogeneous 
reactors.  Using  this  procedure,  an  intermediate 
mechanism  made  up  of  1282  reversible  reactions 
among  282  species  was  obtained.  To  reduce  the 
mechanism  further,  the  flame  structure  was  calcu¬ 
lated  at  one  condition,  and  the  structure  of  the  re¬ 
active  flow  at  two  conditions  close  to  autoignition. 
These  calculations  were  made  using  the  intermedi¬ 
ate  mechanism.  For  all  species,  the  rates  of  produc¬ 
tion  and  consumption  were  computed.  The  results 
were  used  to  determine  the  important  chemical 
paths  of  fuel  breakdown  and  oxidation.  Species  that 
appear  only  along  paths  where  the  rates  of  fuel 
breakdown  and  oxidation  are  small  were  deleted 
from  the  intermediate  mechanism.  This  gave  a  short 
mechanism  made  up  of  770  reversible  elementary 
reactions  among  159  species.  The  intermediate 
mechanism  and  the  short  mechanism  are  available 
in  Ref.  [26], 

Figure  2  shows  the  ignition  delay  time,  rig,  in  stoi¬ 
chiometric  mixtures  of  n -heptane-vapor  and  air  as  a 
function  of  the  reciprocal  of  the  initial  temperature, 
Tu,  for  various  values  of  pressure.  All  numerical  re¬ 
sults  agree  well  with  experimental  data  obtained  in 
shock  tubes  [15]  and  in  rapid  compression  machines 
[16],  Values  of  rig  calculated  using  the  short  mech¬ 
anism  agree  well  with  those  calculated  using  the  de¬ 
tailed  mechanism.  The  agreement  between  values  of 
rig  calculated  using  the  intermediate  mechanism  and 
the  detailed  mechanism  is  even  better.  A  similar  re¬ 
sult  was  found  when  concentrations  of  key  inter¬ 
mediate  species  including  soot  precursors  were  com¬ 
pared  among  the  three  mechanisms.  The  short  and 
the  intermediate  mechanism  were  used  to  calculate 
the  structure  of  the  reactive  flow-field  in  the  coun¬ 
terflow  configuration  and  critical  conditions  of  ex¬ 
tinction  and  autoignition. 

Extinction  and  Autoignition 

Figure  3  illustrates  the  procedure  employed  to  de¬ 
duce  the  strain  rate  at  autoignition  at  a  given  value 
of  T2.  The  maximum  temperature  in  the  flow-field, 
Tmax.  is  plotted  in  Fig.  5  as  a  function  of  the  strain 
rate.  This  figure  is  the  so-called  C-shaped  curve.  For 
values  of  strain  rate  greater  than  <z2,i.  two  solutions 
obtained  are  shown  in  Fig.  3.  The  solution  with  the 
lower  value  of  the  maximum  temperature  is  stable, 
while  the  solution  with  the  higher  value  of  the  max¬ 
imum  temperature  is  unstable.  The  quantity  «2,i  rep¬ 
resents  the  strain  rate  at  autoignition.  A  similar  pro¬ 
cedure  is  employed  to  obtain  the  critical  conditions 
of  extinction. 
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Fig.  4.  Mass  fraction  of  oxygen  in  the  oxidizer  stream  at 
extinction,  Y02>2,  as  a  function  of  the  strain  rate,  a2>e.  The 
symbols  represent  measurements  [12].  The  solid  line  rep¬ 
resents  results  of  numerical  calculations  performed  using 
the  short  chemical-kinetic  mechanism.  The  dashed  line 
represents  results  of  numerical  calculations  performed  us¬ 
ing  the  intermediate  chemical-kinetic  mechanism. 


Fig.  5.  Oxidizer  temperature  at  autoignition,  T2I,  as  a 
function  of  the  strain  rate,  a2, i-  The  triangles  represent 
measurements  [13,14],  The  solid  line  represents  results  of 
numerical  calculations  performed  using  the  short  chemi¬ 
cal-kinetic  mechanism.  The  dashed  line  represents  results 
of  numerical  calculations  performed  using  the  intermedi¬ 
ate  chemical-kinetic  mechanism.  The  plus  symbols  indicate 
conditions  used  in  numerical  investigation. 

Figure  4  shows  the  mass  fraction  of  oxygen  in  the 
oxidizer  stream  at  extinction,  Yo2,2,  as  a  function  of 
the  strain  rate,  a2  e.  The  symbols  represent  measure¬ 
ments  reproduced  from  Ref.  [12],  Fig.  4  shows  that 
the  results  of  numerical  calculations  obtained  using 
the  short  mechanism  agree  well  with  those  calcu¬ 
lated  using  the  intermediate  mechanism.  Numerical 


results  agree  well  with  experimental  data.  The  agree¬ 
ment  between  the  experimental  data  and  numerical 
results  shown  in  Fig.  4  is  significantly  better  than 
that  shown  in  Ref.  [12], 

Figure  5  shows  the  oxidizer  temperature  at  au¬ 
toignition,  T2  j,  as  a  function  of  the  strain  rate,  a2J. 
The  symbols  represent  measurements  reproduced 
from  Refs.  [13,14],  Fig.  5  shows  that  the  results  of 
numerical  calculations  performed  using  the  short 
mechanism  agree  well  with  those  calculated  using 
the  intermediate  mechanism.  Numerical  results 
agree  well  with  experimental  data  at  low  values  of 
the  strain  rate.  At  high  values  of  the  strain  rate,  the 
numerical  results  show  autoignition  to  take  place  at 
higher  values  of  T2 j  in  comparison  to  those  mea¬ 
sured.  The  plus  symbols  in  Fig.  5  indicate  the  con¬ 
ditions  used  in  numerical  investigation  described  in 
the  following  section. 


Influence  of  Strain  on  Autoignition 

Analysis  was  carried  out  to  investigate  the  role  of 
strain  on  autoignition.  Numerical  calculations  were 
performed  at  three  conditions  identified  here  as  case 
A,  case  R,  and  case  C.  Case  A  refers  to  calculations 
at  a2  =  200  s-1  and  T2  =  1164  K,  case  R  at  a2  = 
500  s-1  and  T2  =  1252  K,  and  case  C  at  a2  =  500 
s-1  and  T2  =  1164  K.  Calculations  were  performed 
using  the  short  mechanism  for  YF1  =  0.387,  T1  = 
378  K,  and  Yo2,2  =  0.233.  At  a2  =  200  s_1,  the 
calculated  value  of  T2  I  =  1165  K,  and  at  a2  =  500 
s-1,  T2i  =  1253  K.  Thus,  case  A  and  case  B  are 
close  to  the  critical  conditions  of  autoignition. 

Figure  6  shows  the  maximum  value  of  the  sensi¬ 
tivity  coefficient  for  various  elementary  reactions 
given  by  (k/YOH)(3YOH/d/c)max,  where  YOH  is  the 
mass  fraction  of  OH.  For  all  three  cases,  the  top  four 
reactions  identified  by  the  sensitivity  analysis  are  the 
same,  with  the  chain-branching  reaction  H  +  02  = 
OH  +  O  being  the  most  important  reaction.  This 
chain-branching  reaction  plays  a  pivotal  role  in  many 
applications  where  high-temperature  chemical  pro¬ 
cesses  are  dominant.  Thus,  the  sensitivity  analysis 
indicates  that  high-temperature  chemical  processes 
control  autoignition  in  the  non-premixed  system  in¬ 
vestigated  here. 

It  is  of  interest  to  focus  on  the  differences  among 
the  three  cases  considered  here.  Fig.  6  shows  the 
value  of  the  sensitivity  coefficient  for  the  reactions 
2-C7H1502  =  2-C7H14OOH-4  and  3-C7H1502  =  3- 
C7H14OOH-5  in  case  A  to  be  much  larger  than  those 
in  cases  B  and  C.  They  represent  the  second  generic 
step  in  the  reaction  sequence  that  leads  to  chain 
branching  under  low-temperature  conditions.  The 
first  generic  step  is  heptyl  radical  addition  to  molec¬ 
ular  oxygen  to  form  C7H1502.  This  compound  isom- 
erizes  (above)  to  C7H14OOH,  followed  by  reaction 
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H  +  02  ->  O  +  OH 
CH3  +  H02  ->  CH30  +  0H 
H02  +  OH  ->  H20  +  02 
H  +  02(+M)  ->  H02  (+M) 
HCO  +  M  ->  H  +  CO  +  M 
HCO  +  02  ->  CO  +  H02 
C2H3  +  02  ->  CH2CH0  +  O 
A-C3H5  +  H02  ->  C3H50  +  OH 

2- C7H1502  >  2-C7H1400H-4 

CO  +  OH  ->  C02  +  H 
CH20  +  OH  ->  H0CH20 

3- C7H1502  ->  3-C7H1400H-5 

CH3  +  H02  ->  CH4  +  02 
H202  (+M)  ->  OH  +  OH  (+M) 
CH20  +  OH  ->  HC0  +  H20 
C2H3  +  02  ->  CH20  +  HC0 
C2H4  +  O  ->  CH2CH0  +  H 
C2H5  (+M)  ->  H  +  C2H4  (+M) 
OH  +  H2  ->  H  +  H20 
C2H4  +  OH  ->  C2H3  +  H20 
C2H4  +  O  ->  CH3  +  HCO 
N-C7H16  +  OH  ->  2-C7H15  +  H20 
2-C7H14QOH-4  +  02  ->  2-C7H1400H-4-02 


Maximum  Sensitivity  Coefficient, 


y0H  dk  Jmax 


Fig.  6.  Sensitivity  coefficient  for  various  elementary  reactions.  The  sensitivities  were  calculated  using  the  short  mech¬ 
anism  for  YK1  =  0.387,  T1  =  378  K,  and  Y0z  2  =  0.233. 


Fig.  7.  Profiles  of  temperature,  T,  and  mass  fractions  of 
C7H16,  02,  2-C7H14OOH-5,  2-C7H15,  OH.  H202,  and 
H02.  These  profiles  are  calculated  using  the  short  mech¬ 
anism  for  a 2  =  200,  T2  =  1164  K,  YF1  =  0.387,  7’,  =  378 
K,  and  Y022  =  0.233  (case  A). 


with  02  that  leads  to  the  production  of  two  Oil  rad¬ 
icals  and  chain-branching.  This  is  the  classical  low- 
temperature  branching  sequence  [10],  Fig.  6  shows 
the  value  of  the  sensitivity  coefficient  for  the  reaction 
2-C7H14OOH-4  +  Oz  =  2-C7H1400H-4-02  in 
case  A  to  be  much  larger  than  those  in  cases  B  and 
C.  This  result  clearly  shows  that  low-temperature 
chemistry  is  important  at  low  strain  rates.  Fig.  6 
shows  that  case  C,  which  has  the  same  value  for  T2 
as  case  A,  does  not  exhibit  low-temperature  chem¬ 
istry.  Thus,  strain  has  a  greater  influence  than  oxi¬ 
dizer  temperature  on  low-temperature  chemistry. 

It  is  of  importance  to  examine  the  influence  of 
various  reactions  on  the  concentration  of  OH  be¬ 
cause  it  is  the  primary  reactive  radical  produced  by 
ignition  reactions.  It  is  also  the  primary  radical  that 
consumes  the  fuel  and  intermediate  products  under 
ignition  conditions.  Fig.  7  shows  profiles  of  various 
species  calculated  using  the  short  mechanism  at  a 
condition  close  to  autoignition  (case  A).  Here,  dis¬ 
tances  are  measured  from  the  fuel  boundary.  The 
profile  of  H02  shows  an  interesting  double  peak. 
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Fig.  8.  Profiles  of  the  mass  fraction  of  OH.  f  OH-  cal¬ 
culated  using  the  short  mechanism  for  YF1  =  0.387,  7']  = 
378  K,  and  Y0^  =  0.233. 

The  first  peak  is  clue  to  H02  production  from  low- 
temperature  reaction  paths  given  by  R  +  02  — >  R02 
— >  QOOH  — >  alkene  +  H02.  Here,  Q  represents 
C„H2)j  species  or  structures.  The  second  HOz  peak 
is  due  to  high-temperature  paths  leading  to  H02 
(e.g.,  HCO  +  02  — >  CO  +  H02).  The  profile  of 
OH  shows  a  small  peak  around  6  mm  followed  by  a 
large  peak  around  9  mm  with  a  “shoulder”  around  8 
mm.  The  shape  of  this  profile  is  controlled  by  con¬ 
vection,  diffusion,  production,  and  consumption  of 
OH.  From  a  detailed  reaction  path  analysis,  the  main 
reactions  that  consume  OH  are  found  to  be 


n-C7H16  +  OH  = 

C7H15  +  Ha 

0  - 

45%. 

(Rl) 

ch2o  +  OH  = 

HCO  +  H, 

2o  - 

-13% 

(R2) 

H02  +  OH 

=  h2o  + 

o2 

-7% 

(R3) 

The  main 

reactions 

that  produce  OH  are 

H  +  Oa 

=  OFI  + 

O  17% 

(R4) 

ch3 

+  ho2 

=  ch3o  + 

OH 

11% 

(R5) 

h2o2 

+  M  = 

OH  +  OH 

+  M 

8% 

(R6) 

c7h14ooh  = 

c7h14o  + 

OH 

14% 

(R7) 

Here,  the  percentages  refer  to  the  contribution  of 
the  specified  reaction  to  production  of  OH  inte¬ 
grated  across  the  flow  field.  Reaction  R1  includes  all 
reactions  between  n-C7H16  and  OH  that  give  differ¬ 
ent  isomers  of  C7H15,  and  reaction  R7  includes  all 
isomers  of  C7H14OOH  and  C7H140.  The  small  peak 
in  the  profile  of  OH  around  6  mm  is  due  to  produc¬ 
tion  of  this  radical  from  low-temperature  reactions 
(e.g.,  reaction  R7  where  a  hydroperoxyheptyl  radical 
produces  a  cyclic  ether  and  an  OH  radical).  This 
small  peak  occurs  at  775  K,  which  as  shown  in  Fig. 
2,  is  in  the  negative  temperature  coefficient  region 


(670-830  K)  of  a  stoichiometric  n -heptane/air  mix¬ 
ture  in  a  homogeneous  reactor  at  a  pressure  of  1  bar. 
At  increasing  distance  from  the  small  peak  of  OH, 
the  value  of  YOH  begins  to  increase  as  the  fuel  is 
nearly  consumed  and  intermediate  species  are 
formed.  Prior  to  this  point,  the  mass  fraction  of  OH 
is  small  because  it  reacts  with  the  fuel  via  reaction 
Rl.  Later,  formaldehyde  and  other  intermediate  hy¬ 
drocarbon  species  such  as  ethene  and  propene  dis¬ 
appear,  which  leads  to  a  further  increase  in  the  value 
of  Y0h  for  similar  reasons.  A  number  of  reactions 
produce  OH,  which  includes  its  production  from  de¬ 
composition  of  H202  via  reaction  R6.  Reaction  path 
analyses  at  different  locations  show  that  this  reaction 
gives  the  “shoulder”  in  the  profile  of  Yqh  around  8 
mm.  Reactions  R4  and  R5  make  a  significant  con¬ 
tribution  to  the  rates  of  production  of  OH  around 
7.5  mm,  where  the  temperature  and  mass  fractions 
of  H.  CH3,  and  H02  are  close  to  their  peak  values. 
Figure  7  shows  a  small  increase  of  fuel  around  5  mm 
caused  by  thermal  diffusion.  The  increase  in  tem¬ 
perature  from  the  fuel  boundary  enhances  diffusion 
of  nitrogen  relative  to  the  fuel. 

Figure  8  shows  profiles  of  YOH  for  the  three  cases 
considered  here.  With  increasing  strain  rate,  the 
small  peak  in  the  profile  of  OH  disappears.  This 
again  indicates  that  strain  has  a  greater  influence 
than  oxidizer  temperature  on  low-temperature 
chemistry. 

Concluding  Remarks 

The  short  chemical-kinetic  mechanism  for  n-hep- 
tane  was  assembled  from  the  detailed  mechanism  of 
Curran  et  al.  [10]  without  introducing  any  empirical 
approximation.  For  non-premixed  systems,  this 
mechanism  remains  to  be  tested  at  pressures  greater 
than  1  bar.  It  is  still  too  large  for  application  in  com¬ 
putational  fluid  dynamic  codes.  To  make  such  appli¬ 
cations  tractable,  it  will  be  useful  to  derive  reduced 
chemical-kinetic  mechanisms  made  up  of  global 
steps  from  this  short  mechanism. 

Acknowledgments 

The  research  at  UCSD  was  supported  by  the  U.S.  Army 
Research  Office  through  grant  no.  DAAD19-99-1-0259, 
under  the  direction  of  David  M.  Mann.  Part  of  this  work 
was  performed  under  the  auspices  of  the  U.S.  Department 
of  Energy  by  the  Lawrence  Livermore  National  Laboratory 
under  contract  no.  W-7405-ENG-48.  The  authors  also  ac¬ 
knowledge  support  of  the  German  Science  Foundation 
(DFG). 

REFERENCES 

1.  Warnatz,  J.,  Proc.  Combust.  Inst.  20:845-856  (1984). 

2.  Chevaliar,  C.,  Louessard,  P„  Muller,  U.  C.,  and  War¬ 
natz,  J.,  “A  Detailed  Low-Temperature  Reaction 


2036 


LAMINAR  DIFFUSION  FLAMES— Counterflow  Flames 


Mechanism  of  n-Heptane  Auto-Ignition,”  in  Proceed¬ 
ings  of  the  2nd  International  Sy  mposium  on  Diagnos¬ 
tics  and  Modelling  of  Combustion  in  Internal  Engines, 
S  93-97,  The  Japan  Society  of  Mechanical  Engineers, 
Tokyo,  Japan,  1990. 

3.  Westbrook,  C.  K.,  Warnatz,  J.,  and  Pitz,  W.  J.,  Proc, 
Combust.  Inst.  22:893-901  (1988). 

4.  Bui,  M.,  and  Seshadri,  K.,  Combust.  Sci.  Technol. 
79:293-310  (1991). 

5.  Lindstedt,  R.  P.,  and  Maurice,  L.  Q.,  Combust.  Sci. 
Technol.  107:317-353  (1995). 

6.  Ranzi,  E.,  Gaffuri,  R,  Faravelli,  T.,  and  Dagaut,  P., 
Combust.  Flame  103:91-106  (1995). 

7.  Maurice,  L.  Q.,  “Detailed  Chemical  Kinetic  Models 
for  Aviation  Fuels,”  Pli.D.  thesis.  University  of  Lon¬ 
don,  1996. 

8.  Bollig,  M.,  Pitsch,  H.,  Hewson,  J.  C.,  and  Seshadri,  K.. 
Proc.  Combust.  Inst.  26:729-737  (1996). 

9.  Held,  H.  J.,  Marchese,  A.  J.,  and  Dryer.  F.  L.,  Com¬ 
bust.  Sci.  Technol.  123:107-146  (1997). 

10.  Curran,  H.  J.,  Gaffuri,  P..  Pitz,  W.  J„  and  Westbrook, 
C.  K„  Combust.  Flame  114:149-177  (1998). 

11.  Lindstedt,  R,  Proc.  Combust.  Inst.  27:269-285  (1998). 

12.  Seiser,  R„  Truett,  L.,  Trees,  D.,  and  Seshadri,  K.,  Proc. 
Combust.  Inst.  27:649-657  (1998). 

13.  Seiser,  R..  Pitsch,  H.,  Seshadri,  K.,  Curran,  H.  J.,  and 
Pitz,  W.  J.,  “Experimental  and  Numerical  Studies  of 
Extinction  and  Autoignition  of  n-Heptane,”  Western 
States  Section  F all  Meeting  of  the  Combustion  Insti¬ 
tute,  University  of  California  at  Irvine,  Irvine,  Califor¬ 
nia,  October  25-26,  1999. 


14.  Seiser,  R.,  Seshadri,  K.,  Piskemik,  E.,  and  Lilian,  A., 
Combust.  Flame.  122:239-349  (2000). 

15.  Ciezki,  H.  K..  and  Adomeit,  G.,  Combust.  Flame 
122:239-349  (2000). 

16.  Minetti,  R.,  Carlier,  M.,  Ribaucour.  M.,  Therssen.  E., 
and  Sochet,  L.  R.,  Combust.  Flame  102:298-309 
(1995). 

17.  Griffiths,  J.  F„  Halford-Maw,  P.  A.,  and  Rose,  D.  J., 
Combust.  Flame  95:291-306  (1993). 

18.  Griffiths,  J.  F„  Halford-Maw,  P.  A.,  and  Mohamed,  C., 
Combust.  Flame  111:327-337  (1997). 

19.  Seshadri,  K.,  and  Williams,  F.  A.,  Int.  J.  Heat  Mass 
Transfer  21(2):251— 253  (1978). 

20.  Peters,  N.,  Prog.  Energy  Combust.  Sci.  10:319-339 
(1984). 

21.  Peters,  N.  and  Rogg,  B.  (eds.).  Reduced  Kinetic  Mech¬ 
anisms  for  Applications  in  Combustion  Systems,  Vol. 
mlS  of  Lecture  Notes  in  Physics,  Springer-Verlag,  Hei¬ 
delberg,  1993,  pp.  1-13. 

22.  Pitsch,  H.,  “Entwicklung  eines  Programmpaketes  zur 
Berechnung  eindimensionaler  Flammen  am  Beispiel 
einer  Gegenstromdiffusionsflamme,”  Masters  thesis, 
RWTH  Aachen,  Germany,  1993. 

23.  Curran,  H.  J.,  Gaffuri,  R,  Pitz,  W.  J.,  and  Westbrook, 
C.  K.,  unpublished  results,  2000. 

24.  Lay,  T„  and  Bozzelli,  J.  W Phys.  Chem,  A  101:9505- 
9510  (1997). 

25.  Knyazev,  V.  D.,  and  Slagle,  I.  R.,  J.  Phys.  Chem.  A 
102:1770-1778  (1998). 

26.  Seiser,  R.,  Pitsch,  H.,  Seshadri,  K.,  Pitz,  W.  J.,  and  Cur¬ 
ran,  H.  J.,  Intermediate  and  Short  Mechanism  for  n- 
Heptane,  http :  /  / www-cms .  llnl.gov/  combustion/ com- 
bustion_home.html. 


COMMENTS 


John  Griffiths,  University  of  Leeds,  UK  Are  you  able  to 
account  for  the  diminishing  importance  of  the  low-tem¬ 
perature  reactions  as  the  strain  rate  is  increased? 

Author’s  Reply.  We  are  able  to  account  for  the  diminish¬ 
ing  importance  of  low-temperature  chemistry  as  the  strain 
is  increased.  In  our  analysis,  the  influence  of  a  given  ele¬ 
mentary  reaction  on  the  value  of  YOH  is  presumed  to  be 
a  measure  of  its  influence  on  critical  conditions  of  autoig¬ 
nition.  Fig.  6  shows  that  at  a  fixed  oxidizer  temperature  of 
1164  K,  the  values  of  sensitivity  coefficients  for  reactions 
that  characterize  low-temperature  chemistry  are  higher  at 
the  lower  strain  rate.  This  clearly  shows  the  diminishing 
importance  of  low-temperature  chemistry  as  the  strain  is 
increased. 

• 

Harsha  Chelliah,  University  of  Virginia,  USA.  You 
showed  excellent  agreement  between  experiments  and 
modeling  for  extinction  conditions,  but  comparisons  were 
poor  for  counterflow  ignition.  Is  it  because  of  the  quasi¬ 
steady  ignition  model  used? 


Author’s  Reply,  consider  the  agreement  between  exper¬ 
imental  measurements  of  oxidizer  temperature  at  autoig¬ 
nition  as  a  function  of  the  strain  rate  and  results  of  nu¬ 
merical  calculation  to  be  good.  Close  to  critical  conditions 
of  autoignition,  in  the  experiments  the  oxidizer  tempera¬ 
ture  was  increased  by  1  K.  At  each  temperature,  steady 
conditions  were  maintained  for  approximately  one  minute. 
The  temperature  was  increased  until  autoignition.  It  is 
therefore  appropriate  to  use  the  quasi-steady  ignition 
model  in  the  computations. 


Philippe  Dagaut,  CNRS/LCSR  Orleans,  France.  Why  is 
the  cool  flame  more  pronounced  in  the  model  than  in  the 
experiments  (ignition  of  n-heptane  versus  1/T  in  homoge¬ 
neous  conditions)  at  3.2  and  13  bar  and  not  at  40  atm?  Is 
it  due  to  the  thermochemistry  or  the  kinetics  of  the  data 
obtained  in  this  regime? 

Author’s  Reply.  We  do  not  consider  the  relatively  minor 
differences  between  simulated  and  experimental  ignition 
delays  to  be  important.  These  simulations  were  performed 
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with  the  assumption  of  homogeneous  mixtures  and  no  heat 
loss.  Spatial  variations  in  temperature  inside  the  experi¬ 
mental  device  and  heat  losses  could  affect  the  agreement 
between  the  model  and  experiment. 

However,  your  comment  raises  a  very  important  ques¬ 
tion:  why  does  the  negative  temperature  coefficient  (NTC) 
behavior  become  less  pronounced  with  increasing  pres¬ 
sure?  This  effect  is  observed  both  in  the  experiments  and 
in  the  simulation  (Fig.  2).  As  pressure  is  increased,  the  R 
"h  ()2  «-»  R02  equilibrium  is  shifted  toward  formation  of 
R02.  At  higher  pressures,  the  temperature  must  be  raised 
higher  than  at  lower  pressures  to  dissociate  R02,  a  process 
which  shuts  down  the  low-temperature  chemistry  and 
marks  the  start  of  the  NTC  region.  Thus,  the  temperature 
at  the  start  of  the  NTC  region  increases  with  pressure  (Fig. 
2).  However,  the  end  of  the  NTC  region  is  caused  by  the 
influence  of  the  dissociation  of  H202:  H202  +  M  =  OH 
+  OH  +  M.  The  reactive  OH  radicals  produced  enhance 
the  autoignition  process.  This  reaction  is  pressure  depen¬ 
dent,  so  that  the  rate  of  dissociation  increases,  with  pres¬ 
sure  becoming  significant  at  lower  temperatures  when  the 
pressure  is  raised.  Thus,  at  higher  pressures,  dissociation 
of  H202  terminates  the  NTC  region  before  the  NTC  be¬ 
havior  can  become  as  pronounced  as  that  seen  at  lower 
pressures  (Fig.  2). 

• 

Hai  Wang,  University  of  Delaware,  USA.  Based  on  your 
study  and  a  series  of  previous  work  on  the  ignition  of  C2- 
C4  fuels  in  the  counterflow  configuration,  it  is  quite  con¬ 
clusive  that  ignition  is  caused  by  high-temperature  chain¬ 
branching  chemistry  under  most  conditions.  What  remains 
uncertain  is  the  role  of  heat  release  at  the  ignition  turning 
point.  Have  you  looked  into  the  coupled  effects  of  chain 
branching  and  heat  release  on  ignition  strain  rate? 

Author’s  Reply.  Numerical  computations  carried  out  in 
our  study  include  the  coupled  effects  of  chain  branching 
and  heat  release.  We  have  not  compared  their  relative  in¬ 
fluence. 


Frederick  Dryer,  Princeton  University,  USA.  You  have 
chosen  this  configuration  to  emphasize  the  effect  of  dif¬ 
fusion  on  autoignition.  It  appears  that  while  chemistry  and 
thermochemistry  sensitivities  or  results  are  considered,  one 
might  expect  sensitivities  to  the  assumption  of  particular 
diffusion  coefficients  used  in  the  model.  Indeed,  the  whole 
point  is  that  some  sensitivities  may  approach  those  of 
chemistry  or  thermochemistry.  Have  you  investigated  these 
sensitivities?  Also,  is  the  major  diffusive  sensitivity  only 
thermal  diffusion? 

Author’s  Reply.  We  have  not  carried  out  detailed  inves¬ 
tigations  of  the  sensitivities  of  molecular  diffusivities  of 
various  species  and  thermal  diffusivity  on  critical  condi¬ 
tions  of  extinction  and  autoignition.  As  pointed  out,  these 
sensitivities  may  be  comparable  to  those  of  key  elementary 
reactions.  The  present  study  emphasizes  the  influence  of 
flow  field  and  chemistry  on  extinction  and  autoignition. 
Even  if  molecular  diffusion  is  found  to  play  a  key  role  in 
autoignition  and  extinction,  it  is  unlikely  to  change  the  key 
findings  in  this  paper. 


Chung  K.  Law,  Princeton  University,  USA.  We  are 
pleased  with  your  finding  that  high-temperature  chemistry 
controls  the  ignition  of  strained  diffusive  systems.  It  agrees 
with  the  suggestion  from  our  previous  study  on  the  coun¬ 
terflow  ignition  of  ethane,  propane,  and  butane  at  pressures 
below  5  atm  [1]  that  “ignition  ...  is  initiated  by  fuel  oxi¬ 
dation  following  the  liigh-temperature  mechanism  of  rad¬ 
ical  chain  branching  and  with  little  contribution  by  low-to- 
intermediate  temperature  chemistry.” 
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